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Abstract: We continue the study of hadronic scattering amplitudes at high energy by system- 
atically including nonlinear effects of finite partonic density in hadronic wave function as well as 
the effects of multiple rescatterings in the scattering process. In this paper we derive expressions 
for a single inclusive gluon production amplitude and multigluon inclusive production amplitudes 
when the rapidities of all observed gluons are not very different. We show that at leading order 
these observables exhibit a semiclassical structure. Beyond the semiclassical result, we find that the 
gluon emission has some characteristic features different from the JIMWLK and KLWMIJ limits in 
that the gluons are not emitted independently in rapidity space, but have a correlated component 
with correlation length (in rapidity space) of order one. We demonstrate the consistency between 
this feature of the multigluon observables and the Hamiltonian of the QCD Reggeon Field Theory 
(Hrft) derived in the companion paper We also show that the evolution of these observables 
with total rapidity of the process is generated by Hrft of Q . We discuss whether this evolution is 
equivalent to evolution with Hjimwlk a s far as this set of observables is concerned. 



Contents 



[l]. Introduction [l] 

|2]. High Energy Evolution as Seen From ffl. |3| 

A\ Generalities |3| 



2.2 Definitions 



|3|. The amplitudes Q and the scattering amplitude. [8] 

Inclusive gluon production. [ll] 

|4.1| Single gluon inclusive production. [ll] 

4.2 Double and multi gluon inclusive amplitudes. [14| 



4.3| Evolution with rapidity |1§| 
[5]. Discussion. 2C 
|A]. Appendix - Single gluon inclusive amplitude. 23 



1. Introduction 

This is the second of two papers devoted to the study of hadronic observables at high energy taking 
into account the effects of Pomeron loops. In the first paper Q, we derived the Hamiltonian of 
the QCD Reggeon Field Theory Q which accounts for the effects of nonlinearities in the hadronic 
wave function as well as multiple scattering effects in hadronic scattering, and thus fully includes 
Pomeron loops. It supercedes the JIMWLK |§ , [§ , f§ , [§ and KLWMIJ Q evolutions which are 
only applicable in parts of the parameter space. This Hamiltonian - Hrft governs the evolution of 
many hadronic observables with energy. In particular it directly gives the evolution of the forward 
scattering amplitude. The derivation of was based on the earlier results ^ for the evolution of 
the hadronic wave function. 

In this paper we extend our analysis to other observables. In particular we derive expressions 
for the amplitude of inclusive gluon production at high energy. We consider single gluon production 
as well as multi-gluon production when all the emitted gluons are close to each other in rapidity 
so that the rapidity evolution between them does not have to be considered. There has been a 
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lot of interest in the calculation of inclusive gluon production in the framework of the high energy 
evolution. Starting with which calculated single gluon production for DIS in the dipole model, 
this observable (as well as multi-gluon production) was considered beyond the dipole model limit in 
DIS in 0, fU], 0, 0, 0, |15], [T§ and in nucleus-nucleus collisions in ][T7], |§, |l9|, Some, effects of 
the Pomeron loops which do not extend to the rapidity of the measured gluons in DIS setting were 



studied in [21 



In the present paper we extend these results by including the Pomeron loop effects using the 
methods of Q. Our calculation is applicable to both the DIS situation and the nucleus-nucleus 



scattering. We follow the approach and techniques developed in [13, 14]. The extension of these 
results to include evolution between the observed gluons is straightforward but we leave it for future 
publication. 

We note that very recently inclusive gluon production in nucleus-nucleus collision was discussed 
in [ 20 1 . Our present work has many parallels with |20| even though the methods we use are quite 



different. Our results are more complete in the sense that we provide an explicit form of the ob- 
servables in terms of the color charge densities of the incoming hadrons (nuclei). It is not necessary 
in our approach to solve dynamical equations of motion in order to find the outgoing classical color 
fields, as is required by the procedure of |2(J. All that is necessary is to average a function of the 
incoming charge densities (albeit a complicated one) over the initial probability distribution. The 



element of averaging over the initial distributions is also present in the formalism of [20|. 

We also derive a set of subleading corrections to the double inclusive and multi-gluon ampli- 
tudes, which exhibit short range rapidity correlations between the gluons in the final state. These 



correlations are formally subleading in a s and are therefore not present in the analysis of [20|. In our 
approach these terms are however under control. Since our derivation of the multi-gluon observables 
is closely linked with the derivation of Hrft we are able to show quite generally that the expression 
for Hrft derived in |IJ necessarily requires the presence of such short range gluon correlations in 
the multi-gluon observables. We stress that these are not correlations due to the evolution between 
the rapidity of the gluons [l2| , |14| ], whose perturbative tail we also see in this calculation, but rather 
correlations bunching (or anti-bunching) all emitted gluons inside a rapidity interval of order unity. 

Our approach allows us to show immediately that the evolution of multi-gluon observables with 
respect to "global" rapidity parameters is given by the RFT Hamiltonian Hrft derived in By 
" global" rapidity parameters we mean either the total rapidity of the process at fixed gluon rapidity, 
or the rapidity difference between the emitted gluons and the projectile; or between the emitted 
gluons and the target; but not the rapidity differences between the emitted gluons themselves. 

Finally we discuss the question to which extent this energy evolution of the multi-gluon inclusive 



observables can be approximated by the JIMWLK evolution, as suggested in [2C]. Formally it 



appears that Hrft can be replaced by Hjjmw lk as long as one considers the evolution of this set 
of observables, with corrections to Hjimwlk being suppressed by powers of a s . However we give an 
example of a situation where such a formal argument fails, and thus advise caution on the subject. 
We believe this question merits further investigation. 

We start with recapping the notations and some formulae from Q which we will need for our 
derivations. The next section does not contain any new results and is included purely in an attempt 
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to make the present paper self-contained in terms of notations and definitions. All the formulae 
presented in this section are also found in Q. Reader who is familiar with notations and results of 
H| can skip directly to Section 3. 

2. High Energy Evolution as Seen From 
2.1 Generalities 

We consider a hadronic projectile moving to the right with large energy. We separate the degrees of 
freedom of the hadronic wave function into "valence" - gluons with rapidity greater than some fixed 
value, and "soft" - with rapidities lower than this separation rapidity. The properties of the valence 
component of the wave function are characterized by correlators of the color charge density operator 
j a (x). The wave function of the hadron can be written as 

|*) = Sl[a,<J,j\\v) (2.1) 

where a and a* are soft gluon creation and annihilation operators. The valence state \v) has no soft 
gluons and is therefore annihilated by the soft gluon annihilation operators 

a\v) = 0. 

The evolution operator O is a unitary operator of the Bogoliubov type 

O = CB. (2.2) 
Here C is a coherent operator that creates the "classical" Weiszacker- Williams field 

C = ex P |iv^|^ b *( k ) J ^ + 4 a (^)]} • (2-3) 

and B is a Bogoliubov type operator responsible for the leading quantum corrections. The rapidity 
variable is defined as r] = In and the creation and annihilation operators are canonical in the 
rapidity basis 

[attv,k),af(Z,p)} = (2vr) 3 %-0<5 2 (A ; -p). (2.4) 

The Weiszacker- Williams field W(x) depends only on transverse coordinates x and is a two 
dimensional pure gauge field: 

b? = -- f abc U ]bd diU dc (2.5) 
9 

where f abc are the structure constants of the SU(N) and U ab is an SU(N) group element in the 
adjoint representation. The Weiszacker- Williams field is related to the valence color charge density 
by 

«(x) =f(x) (2.6) 

The soft gluon operators a{rf) are labeled by rapidity variable r\ as well as transverse coordinates (or 
momenta), color and rotational indices. They populate the rapidity interval Y equal to the parameter 
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of the boost transformation that brings the wave function of the hadron from \v) to We do not 
denote this cutoff explicitly in most of our formulae and extend the rapidity integration over r] to 
infinity for all quantities for which the integration converges. The Y dependence is important only 
for divergent quantities and eventually this very dependence determines the evolution of physical 
observables with rapidity. 

The explicit form of B is not important here, but its action on the gluon creation and annihilation 
operators is linear 



P a = Ba a B^ = & a/3 \j] ap + $ a p[j) oj, 0+ = Ba\ & = Q* af3 [j] a\ + KpUH ( 2 - 7 ) 



Pi =BalB^ 

It turns out to be useful to define separately the matrix 

N(p,k) = @(p,k)-<f>(p,-k), 
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where we have indicated explicitly only the transverse momentum dependences. The transformation 
matrices 0, <3? and ./V depend on the valence color charge density and where calculated in Q. The 
explicit form of O and <3? is not important for us in this paper. The expression for N is 



N{p, V ;k,0 



-«-»7) 



ie 



-(p\(l-2l)\k)-(p\(l-2L) 



\k) (2.9) 



where D is the transverse covariant derivative in the background filed b, and the longitudinal pro- 
jectors I and L are defined as 

ed 



ah 



djdj cab, 

Q2 > 



L? 3 EE DT (b) 



1 



D 2 (b) 



Df{b) 



The vacuum of the operator (3 as calculated in M is 



\0)p = B\0) 



with 



e 4 



A 



i 1> hi(l - At A) -lotAot 



10) 



(2.10) 



(2.11) 



(2.12) 



Eqs.(2.1. 2~3l , 2.11 ) characterize the wave function of a hadron with rapidity Y. This hadron, 
which we call the "projectile" scatters on another hadron- the "target". The target is specified 
by a distribution of color fields a a . Calculation of many interesting observables in the scattering 
process, including multi-gluon inclusive amplitudes, can be represented as the double average over 
the projectile and the target wave functions of some operator O 



(O)t 



DaW T [a] O ; 



6 = (v\6[j]\v) 



1 



DjW p \j]0\j}. 



(2.13) 



Here is a weight functionals representing the probability distributions of the projectile color 

charge j, while VF T [a] is the same for the target fields a. The rapidity evolution of this observable 
is given by Hrft ( 2.32 ) in the following sense ||], [jj 



dY 



-.0 



Dj W p [j] H RFT 



3, 



Sj 



0[j] 



(2.14) 
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or equivalently 



dY 



Hrft 



(2.15) 



The Hamiltonian Hrft is formally defined as the expectation value of the unitary operator R a in 
the hadronic wave function Iv?) 



Hrft 



Sj 



d 

dY 



(0|Qt[j, a, at] R a n[j, a,<J]\0)\ Y=0 ■ 



with 



and 



J^ a = g fx ■5 , so/t( X ) Sj a (x) 

fsoftW =9 f ^a^ b (rj,x)T b %a c ( V ,x) 



(2.16) 



(2.17) 



(2.18) 



The operator R a shifts the valence color charge density by the charge density of the soft gluons. This 
reflects the fact that in the boosted wave function not only valence, but also soft gluons participate 
in the scattering, and thus the effective color charge density must include both j and j so ft- 

Hrft wa s calculated in [ffl]. To summarize the results of (lj we have to introduce (alas - at the 
risk of being repetitive) some additional notations. 

2.2 Definitions 

First, since the color charge density operators do not commute, we need to keep track of the ordering 
in eq. (|2.16j ). As a result we have to define the right and left color charge density operators Jr and 
Jl. The labels refer to their order in the calculation of the matrix element of O and the Hamiltonian 
Hrft- The operator fi in eq.( |2.16| ) has to be understood as a function of the right charge operator 
Jr, while the operator fit is a function of Jl- This is important, since observable O generically 
depends on the color charge density, and in the calculation of its matrix element in |^) it always 
appears sandwiched between fit and fi. Defining Jr and Jl specifies the relative ordering between 
the color charge operators in the wave function and those in the observable. As shown in |22|| , the 
right and left charge densities are related by the rotation with unitary matrix R (here t is the ordering 
variable @) 



R" 



V exp{o / dt T 



5j c {x,t) 



ab 



(2.19) 



The matrix R has the meaning of a scattering matrix of the gluon of the target on the projectile 1 22]. 



As discussed, for example in [23], it can be taken as the basic degree of freedom of the Reggeon Field 



Theory. We can then write [22] 



gJ a dx) 



-gti ^R(x)T 
-gir \ T a R(x) 



&Rt(x) 
6 

dRT(x) 



(2.20) 



[R(x)J R (x)] a . 
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Accordingly we define two classical fields, 6r and bi: 

b% = -- f abc U^ bd [J R } di U dc [jR] = -- f ahc U ] R d d t C/f ; 

b% = -- f abc U^ bd [J L ] d t U dc [J L ] = -- /^tff* $ U dc (2.21) 

The matrices 0, <3? and A also "bifurcate" into Right and Left quantities. 

A R = A[J R ] = Q-^Jr] <S>[Jr}; Al = A[J L ] = Q'^Jl] $[Jl] • (2.22) 

It is thus clear that for the purpose of the calculation of the average of any observable the proper 
ordering of the factors of the charge density is equivalent to the substitution 

nt_> n t =B t c t. n^n R = c R B R (2.23) 

where the subscript L (R) indicates that the respective operator depends on Ji (Jr). Due to the 
presence of the operator R a in eq. (|2.16D (and in other similar averages) it turns out that the Left 
quantities are always rotated by the matrix R. It is thus convenient to introduce the barred quantities 

n L = C L B L ; B L = R\B L , C L = R\C L R a . (2.24) 

and accordingly 

b a Li = R ]ab b b Li - A L = R^A L R; N L = N L R (2.25) 
Finally we introduce another set of matrices 

^ = (eet-$$V E afi = (d>e T - e3> T ) Q/3; (2.26) 

3 = K~ X E. (2.27) 

which satisfy 

KE T - EK T = 0, KK ] - EE ] = 1; K ] K - E T E* = 1; K ] E - E T K* = ; (2.28) 

1 K*K] = K T K*. (2.29) 



Above we have omitted the indices L and R. All the barred quantities depend on Jl while the 
unbarred ones - on Jr. We will stick to this practice in the following, whenever it does not lead to 
confusion. 

The matrices K, E and 3 are the analogs of the Bogoliubov transform eq.(2.7) for transformation 



between the operators (3r ( defined by eq.( |2.7D with Or, etc.) and (5l ( defined by eq.( |2.7D with Qr, 
etc.). The matrices K and E, which we will need in the following have been calculated in |i 

K(p, V ;k,0 = i(p\2(l - I - L R )R\l -L L ) + (l-l- L R )(A - A R )R<<(1 - I - L L )\k) 

E(p, V ;-k,0 = i(p|2(l - I - L R )R\l -L L ) + (l-l- L R ){A - A R )R\l - I - L L )\k) 

+ t , " 22)^(1 " 2Z)|fc> - i{p\{l - 2L R ) i i?t(l - 2L L )\k) (2.30) 
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with 



D = R ] D L R; 



A = R^AlR. 



(2.31) 



In eq. ( |2.30 ) A is an operator in the transverse space, which we do not specify, since it does not 
appear in the following, D R and Dl are transverse covariant derivatives in the background fields b R 
and &l respectively, and longitudinal projectors are defined as in eq. (|2~10|) . In eq.(p0|) all D R 

are ordered to the left of Dl. 

2.3 Hrft 

The result of O for H R ft to leading order in a s in the concise matrix notation is 
d ( 



Hrft 



dY 



[ [b(x)-b{x)]f N\7 1 ,x;ty)K- 1 (Z,y;\,z)N(\,z;C,u) [b(u)-b(u)] 

Jx,y,z,u Jri^^XC, 

(2.32) 



1 



+ -Tr In (KK*) \\ Y=0 



The Tr In term in this expression is subleading in a s and will not be considered any further. While 
b is independent of rapidity, both N and K depend on rapidity differences only. Introducing 



N ± (x,y) = / N(x,r),y,S) = [l — l — L] 



(2.33) 



which will reappear in the next section and performing the rapidity integrations the result for the 
first term in eq.(2.32) can be recast in the form 



1 



H RFT = ~[b R R ] -b L ) (l-l-L L ) 

IT 

x (1 - 21) R) (1 - 21) + (1 - 2L R ) fit (i _ 2L L j\ ^ (I - I - L R ) [b R - b L ] (2.34) 
Reinstating explicitly transverse coordinate dependences this reads 



H RFT = ^3 / [b b m (x)R^{x) - b a Li (x)] 



x,y,z,z 



1 



[Ul(x)U L (z)} ac K^ k cd (z,z) 



13 (x — z) 2 

1 



Ski 



(x - z)j (x - Z)j 
(x - z) 4 

, (y - z)k (y - z)i 



(y - z) 2 

^ + [U R (z)U R (y)r] [b Rl (y) - R^(y)b{ k {y)) 



(y - z) 4 



(2.35) 



with 



K ± %(x,y) 



2tt 2 



(x - zy 



(x - z)j (x - z) k 
(x - z) A 



x \ eP*{z) + \u R (x)U R (z)R\z)uUz)U L (y) 



ab 



(z - y) k (z - y)j 



{z-yf 



(z-y) A 



(2.36) 



With all the definitions in place we can now proceed to calculation of the multi-gluon spectrum. 
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3. The amplitudes Q and the scattering amplitude. 



We will start by rederiving the results of ]|] in a slightly more general way which later will allow us to 
calculate additional observables. First, we remind the reader the general formalism of p3|, which we 



will use here. The matrix element eq.(2.16) pertinent to the calculation of Hrft can be represented, 
using the optical theorem as 



(o|i-nt£ a o|o) = I(o|(i-fit£tn) (i-rtR a n)\o). 



(3.1) 



This assumes that the amplitude is real, which in our case can be verified directly by examining 
eq.(2.34). Introducing the complete basis of intermediate states and defining the amplitudes 

Qn(xi,Vi) = (xi,TH',...;x n ,ri n \(l-rtRAn)\0) = (n\&n\0) (3.2) 

we can write 



Hrft 



1 d 

2 dY 



71=1 



£ / n™ =1 [dxidm] QU 



(3.3) 



As explained in , the amplitudes Q n has the meaning of emission amplitudes of n gluons in one 
step of evolution. The knowledge of Q n in principle allows one to construct also exclusive observables, 
for example constraining the number of particles produced in the final state. In the next section we 
will use Q n to calculate single and double inclusive gluon production. 

Our aim is to write the expression for Q n in the normal ordered form, so that the amplitudes 
can be easily read off. 



fitniO) = 8*tfCB\Q) = B f e^b-fHa+at) BB*B\0). 



where 



V2(b - b) (a + a f ) = V2 / d 2 x 



b«(x)-b?(x) 



drj 
2^ 



a?(x,r))+af(x,ri) 



(3.4) 



(3.5) 



For the first factor, using the fact that B is the Bogoliubov operator with the action similar to eq.(|2_. 
with 0, $ — > 0, $ we have 

g\ e iV2{b-b) (o+ot) g _ e iV2{b-b)(m a+N T a<) _ ^— (6—6) iVt f} (b—b) g*\/2 (6—6) N T at g*V^2 (6—6) N^a g\ 



The state B^ B |0) is the vacuum of the operator 



= B^BaBB^ = Ka + E, 



(3.7) 



with K and E defined in eq.( |2.26 ). Comparing eq. with eq. (2/7), we can use eq. ( 2.11] ) with 

the substitution A —* S to write 



i Tr ln(l-HtH) 



e 2 



|0) 



Now using the fact that the annihilation operator can be represented as a 



<5at 



we can write 



jV2(b-b) Aft a 



-i n .tR„t 



e 2 



|0> 



-\ (at +iy^(6-6)JVt) H (at +i-\/2 (6—6) Aft ) | q\ 



(3i 



(3.9) 
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And so, finally the normal ordered form is 

fitfilO) = e sTrhi(l-HtH) e -(b-b)[NW-mSN*](b-b) £ iV2 (6-6) [N T -Nt s] at e -i tS t|g\ (3 1Q ) 



Using the relations eqs. ([2^2S| , ^27| ) and the definit ion of H in terms of G, <I>, 0, <I? it is straightforward 
to show 

N = KN - EN*; N = K ] 'N + E 1 ' N* . (3.11) 

Thus we can rewrite eq.( |3.10| ) as 

O f O|0) = e-l^HK^K) e - (6-6) JVt K- 1 N (6-6) g i ^2at A-i JV (6-6) e -£ a t S at| ^ (3^) 

Taking the overlap of eq. fl3.12|) with the vacuum ((0| 0^0 |0}) we reproduce the result of Q eq.( |2.32|) , 
while the last two factors in eq. (|3.12|) allow us to calculate Q n in a straightforward manner. 

Scattering amplitude once again. 

Before moving on to discussion of inclusive gluon spectrum, we rederive the expression for the scat- 
tering amplitude and for the RFT Hamiltonian using eq.( |3,3j ). This will explicitly establish the 
equivalence of the two approaches and will also serve as a consistency check on our expressions for 
the amplitudes Q n . According to eq.(|3.3j) we need to extract the term linear in Y from the sum 
of squares of the amplitudes Q n integrated over the transverse coordinates and rapidities. The am- 



plitudes Q n are obtained by expanding the exponent in eq.(3.12) to appropriate order. In general 
therefore the calculation will involve expressions of the type 

<0| [(6 - b) K^ 1 a] m [a E* a] k [a) K' 1 N (6 - b)} 1 [a) E a f ] n |0> . (3.13) 

Here m, n, I, k are some integer powers to be summed over. For the purpose of the present discussion 
we will think of the matrix H(£ — 77) as short range in rapidity space. Although this is not really the 
case, as discussed in detail in jjj], the nonlocal in rapidity terms have to be subtracted from the final 
result, since they lead to higher powers of Y. Recall that the field b — b does not depend on rapidity. 
Thus after contracting all a's with a^'s in eq.( |3.13| ) we will generate expressions of two types and 
their products (for simplicity of notation we drop in this schematic discussion the factors of N K^ 1 
which accompany powers of (b — b)) 

(6-6) / pE] k (v,Z)(b-b); Trp E} 1 . (3.14) 

It is easy to see that each one of these terms is of order Y . Thus we only need to consider contractions 
m eq. ( fTp that lead to appearance of only one such "connected" term, rather than product of two 
or more - the "disconnected" graphs . This has an immediate consequence that only terms of order 
(6 — b) 2 have to be kept. Higher powers of b necessarily lead to disconnected graphs and thus 
necessarily to higher powers of Y. On the other hand, the terms which do not contain any powers 
of b are suppressed by a power of coupling constant a s . Thus we are lead to the conclusion that the 
only terms that contribute to Hrft are connected contractions of order (b — b) 2 . With this in mind 
the calculation becomes straightforward. 
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Let us consider odd and even n's separately. 



Q 2n+1 = i (-1)" J^L (n| [at A^ 1 N (b — 6)] [«t 3 at]" |0> . 



(3.15) 



Then 



dx i dr 1i Ql n+1 (x,rj)Q 2n+1 (x, V ) = 2 (6 - b) K^ 1 [E ~tf K' 1 N (b - b) . (3.16) 



Thus 



QL+i Q2n+i = 2 (6 - b) iVt ift-i _ _ T K-i at (6 - 6) = 2 (6 - 6) N (6 - b) (3.17) 



n=0 



where we have used eq.( |2.29| ). For even n we have 



Q2n = (-l) n 1 ^ 1 (n\2n[a^K- l N(b-b)] 2 [^Ea^ n - 1 + [a^Ea^] n \0). (3.18) 



2«n! 



So 



dx i dr H Ql n (x,r l )Q 2n (x,r l ) = (b - b) N T K~ 1T ~) [E E^ 1 K~ x N {b - b) + h.c. . (3.19) 



Thus 

oo 



2d 



n=l 



b) N T R- 1T S f - — _^ if- 1 i\T (6 - fe) + (i-^r't - __ 5ji:-i«jv« (6-6) 



(6 -6) N T K~ 1T N(b -b) + (b — b)N^E K~ u N* (6-6) 



(3.20) 



Altogether 

oo 

YjQnQn = (6-6) bjV+JV + N T R- 1T E ] N + A^AT" 1 * aH (6-6). (3.21) 

n=l 

Finally using eq.( |3.11| ) it can be recast in the form 

oo 

Qn Qn = (b - b) [#t K~ l N + /i.e.] (6-6) = 2 (6 - 6) K' 1 N (6 - 6) (3.22) 

n=l 

which reproduces Hrft - the first term in eq.( 2.32| ). 

We note without proof, that were we tc 
also reproduce the second term in eq.( [2.32 ). 



We note without proof, that were we to keep also subleading in a s terms in Q 2n Q 2n we would 
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4. Inclusive gluon production. 



4.1 Single gluon inclusive production. 

We now turn to derivation of the amplitude dn(r],k)/dr] for inclusively produce gluon with rapidity 
r\ and transverse momentum k. Here we follow the approach of (13, 14]. 



The single gluon inclusive amplitude is defined as 

dn(r], k) 



J DjDSW p [j]O g (k, V )W T [S] (4.1) 



dry 

with the gluon operator to be measured in the collision 

O g (k,ri) = — (0| fit St fi a\ a (rj, k) af(rj, k) fitSfi|0) = 

= — / e *H*-v) (Oltf&n a\ a (ri,x)af(r),y) fitSfi|0). (4.2) 

2^ J x ,y 

The state |0) as before is the vacuum of the soft gluon Hilbert space. The operator S here is the second 
quantized S-matrix operator, which acts as a color rotation on the gluon creation and annihilation 
operators as well as on the color charge density in the operator fit 

St a a {x) S = S ab (x) a b {x); S f f(x) S = S ab (x) j b (x) . (4.3) 

The unitary matrix S(x) is the scattering matrix of a single gluon in the color field of the target. 
The target average in eq.(fO^) amounts to integrating over S(x) with a weight function determined 



by the wave function of the target E^],] 14]. 



The definition eq.(4.2) involves objects almost identical to the amplitudes Q n discussed above. 
In fact the operator fit fi becomes identical to the operator fitSfi upon substitution R — > S. The 
operator O g is computed in the Appendix by acting with fi directly on fields a. Here we will use 
a different procedure which makes clear the connection to the calculation of Hrft in the previous 
section. Inserting resolution of identity in eq. (|4.2| ) we can write 



°<?(M) = ^fjj[ z ik(x ~ y) J j^^I^-l d 2 z l d Vl nQi[S-,z l ,r Jt] x, V }Q n [S-,z uVl -,y,r ] ] (4.4) 

where Q n [S] is obtained from the amplitude Q n by substituting S for R. In the rest of this section 
we will simply use Q n to denote the amplitudes that enters eq.([4.4j), but will keep in mind that they 
are related to the amplitudes used in the previous subsection by the aforementioned substitution. 
Note that now Q n are operators on the target Hilbert space and thus the averaging over the target 
as in eq.(4.1) is necessary. 



The expression in eq.( |4.4j ) can be calculated using eqs,( |3.16 , |3.19|) . In practical terms we have 



to cut one index linking one of the S's and E)t (or S and b — b etc.) in eq. ( jTilj ) and eq.(|I|) and 



take the free transverse coordinate of H equal to x, while that of Ht equal to y. Each term can be 
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cut in n possible positions, and all those cuts have to be summed over. Thus for odd n we obtain 
J2 (2n + l)Ql n+1 (z 

ii 

,„_q J dzi,rn 



n=0 
n-l 



^2 (b-b) N[ Kt- 1 [E ~X *[ S St ]"" m K ~ l N ± (b-b) + 



m=0 



+ £ (5 - b) N{ A-t-i [h H t r ~ x yH t [3 Si 



in— m— 1 



JV ± (6 - b) 



m=0 



= 2(b-b)N{ rt-^l-ESitf ^l-EE^K- 1 N ± (b-b) + 

+ 2(b-b)N[ ift- 1 [1 - S St]- 1 S x ,3t [1 - S St]- 1 K- 1 N± (b-b) = 

= 2 (b-b) N[K y X K^ N ± (b-b) + (b- b) N[ E y x Eft N± (p _ b ) (4.5) 

where we have used a shorthand notation x ^46 = J A(x,y)b(y); bA y = j x b(x)A(x,y). 
For even n we obtain 



00 . 

^2 2n Ql n (zi,Vi, x ,v)Q2n(zi,rji,y,rj) 

n=l ^ z "^ 



£ 

n=\ 
n-l 



J2(b-b) Nl K T - X St [S St]™ x [S St]™"™ K- 1 N± (b-b) + 



m=0 



+ (b - b) Nl K T ~ X [St S] m ,St [S S t r "- 1 K -i N± { l-b) + h.c.(x - y) 



m=0 



(S-fcjivlJ^-Mi-sTSl-, 1 ^!! 



tm-i ~t 



1-1 6—1 



R-'Nx^-b) + 



+ (b-^iVj^-^l-Sts]- 1 ^Stfl-SSt]- 1 ^- 1 ^^-^) + h.c.(x^y) 
= 2(1 - b) Nl E* X K ] N ± (b-b) + 2 (b - b) n[ E y x K T N* ± (b - b) 



(4.6) 



Adding eqs.( f4.5p and ( |4.6| ) and using eq.( p.ll ) it is now straightforward to show that 



7T 



7T 



O g [k, V ] = - I e lk ^(b-b)m x N ± (b-b) 



x,y 



x.y 



J- v 



J k (x-y) 



(b L - b R 5t) (I -I- L L ) y x (l-l-L L ) (b L -Sb R ). (4.7) 



Note that the gluon emission operator O g is independent of the rapidity of the gluon 77, as it should 
be for a boost invariant plateau. 



As a corollary we note that integrating eq.(4.7) over the transverse momentum we obtain the 
expression for total multiplicity 



dn 1 
dij it 



( (b L - b R St) [1 - l - L L + L L l + IL L ) (b L - S b R ) ) j>3 . 



(4.8) 
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This generalizes the results of j|] for the inclusive gluon spectrum in p-A scattering. We can check 
explicitly that in the KLWMIJ limit eq. (|4.7| ) indeed reduces to the known result of (9|, [R], [n], |l3| . In 
the limit of dilute projectile as before we take 



L^l; ^(x)^i/^|f^5-/(,). (4.9) 



In this limit then the two factors 1 — 21 = 5 l i — cancel against each other. We then obtain 

= A / J k{x - y) ( j\z) \s a \z) - S ab (y)} [S ac (z) - S ac {x)\ f{z) ) hS 

"'V 47T Jx,y,z,z L - 1 \U z ) \ x z ) 

(4.10) 

which is the known result [ll], 13]. 



Recently single inclusive gluon spectrum in nucleus- nucleus collision was discussed in [ pC|| . Al- 
though the setup of [pO] is somewhat different from ours, it is possible to establish close correspon- 



dence between the two approaches. The procedure of [20] is the following. The single gluon inclusive 
spectrum is defined as 



dr\ 7r 



)p,t (4.11) 



where A is the solution of classical Yang-Mills equations of motion with initial condition corresponding 
to colliding sheets of color charge density (nuclei) |24| . The solution should be taken at asymptotically 
large time after the collision. The averaging is then done over initial conditions with separate weight 
functions for the projectile and target color charge distributions. This procedure is somewhat implicit 
since the classical field A has to be found by numerical solution of the classical equations. 



Comparing this to our result, we see that eq.( |4.7| ) indeed '"feels'" like the classical field. The 
expression 

A = N ± (b - b) = (1 - I - L L ) [b L -Sb R }. (4.12) 

plays the role of the classical field at asymptotically large times after the collision. With this iden- 
tification our formula becomes very similar to that of |2(| except for the fact that our expression is 
rather more explicit. To evaluate eq. ([4.7| ) there is no need to solve dynamical equations of motions, 
but rather only to solve the static classical equations at early time, which express the classical field 
b in terms of the color charge density j. In fact it may be possible to avoid this step altogether if 
one can write down directly a weight function for b as advocated in the fourth paper in Q. This is 
an interesting question and we intend to come back to it in the future. 



We also mention, that an expression for dn ^t k ^ was suggested in [17]. The averaging over the 



target and projectile weight functionals was performed in [17] using the McLerran-Venugopalan model 



1 25]. Since we are still not at the point where we can calculate the average of the expression (|4.7| ), 
we cannot compare our result with the suggestion of [Oj. 



Ref.[20| also discusses the evolution of the single gluon spectrum with energy. We have so far 
discussed only the explicit form of the observable itself. The evolution in our approach however is 
very straightforward to understand. We will discuss this later on in this section. 
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4.2 Double and multi gluon inclusive amplitudes. 



The same strategy can be used to calculate the double inclusive cross section as well as higher gluon 
number correlations. In this paper we only consider multi gluon observables where all the counted 
gluons have rapidities not too far away from each other, so that it is not necessary to consider rapidity 
evolution between them. The double gluon cross section for two gluons with rapidities rj and £ and 
transverse momenta k and p is defined as 

dn(V d ^ ,P) = J^((0\^ S^4 a ( V ,k)aUri,k)4 a ^p)aU^p)^ Sn\0)) hS (4.13) 



(27 



^ / e ik(x-*)+i P (y-y) ( (0 | n t £t n t« {rj} x) a a {r]} 2) fl t« ^ y) a » {U) rtsn\ 0) 

J(x,x;y,y) 
I ■ 

^2 



Jk(x-x) + ip(y-y) 



n=2 x > x > 



— U^d 2 z i dr ]l n(n-l) 



(Qn(S;Zi,rn;x,r],y,C) Q n (S;zi,r)f,x,r),y,£) ) jjS ■ 



The calculation is straightforward albeit fairly long. The leading contribution to the double inclusive 
amplitude is obviously O(^). The result for this contribution is 



dn(rj,k;£,p) 
drjd£ 



1 

71-2 J(x,x;y,y) 



Jk(x-x)+ip(y-y) 



(( (6 - b)N{ s x N ± (b -b) (6 - b)N{- y N ± (b - b) ) hS 



= \\ e *l*-z)+rtv-v){\b L - b R tf)(l - Z - L L ) S x {l-l- L L )(b L - Sb R ) 

11 J{x,x;y,y) 

~b L - b R tf)(l - I - L L ) y y (l-l- L L )(b L - Sb R )] ) jtS ■ 

(4.14) 

This again is the "classical" contribution. It is "disconnected" in the sense that it does not depend 
on either rapidity rj or £ and as an operator on the valence and target Hilbert spaces, is equal to 
the square of the single inclusive amplitude O 2 . This result immediately generalizes to multi-gluon 
inclusive amplitudes defined as 

dn(r]i,ki; ...,7] n ,k v 



dr)i...dq n 



= j DjDSW p [j]T n (k 1 ,...k n ; Vl ,... Vn )W T [S} (4.15) 



with 



r n (fci, ...k n ;r]i, ...rj n ) = (0| & SI aj a (^, h) a?( Vl , h) ...af(r, n , k n ) a^, k n ) SSl\Q) . 

(4.16) 

Similarly to the double gluon case, T n can be expressed in terms of Q n : 



J ki (xi—xi) +...+ik„ (x 



-x ) f 1 n! 

J (2vr) m - n ml 



1 f 

x IT™~ n d 2 Zi drji Qj n (S , ;z i ,^;si,»7i,...,x n ,%) Q m (S; z h rjf, x u rji, ...,x n ,rj n ) . 
To leading order in the coupling constant 



Y n {k 1 ,...,k n ;r} 1 ,...r] n ) -> T n {k 1 ,...,k n ) = O g (k 1 )...O g (k n ) 



(4.17) 
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or equivalently 

dn , (fcl, 7 few) = - n [ elklixi - Sl)+ ... +i k n ( Xn - Sn){[A{xi) . A( x l)UA{Xn) . A{2n)] ks 

ar\\...ar\ n tt J ( X i,x 1 -,...x n ,x n ) 

(4.18) 



with the "classical field" A defined in eq.(4.12). Thus all these observables have the same structure 
as discussed in p0| with the identification of the classical field at asymptoticaly late times given by 
eq.( f4.12| ). All these leading order contributions do not depend on the rapidities of emitted gluons, 
since the gluons in leading order are emitted independently. This is the same as in the leading order 
BFKL calculations [^] . This of course does not mean that the multi-gluon amplitude is simply the 
product of single gluon ones 

dn{ki, ...,k n ) ^ dn(k{) dn{k n ) 
dr)\...dr] n <% dr) n 

It is the averaging over the projectile valence and the target Hilbert spaces that breaks the factor- 
ization, even though the equality holds on the level of operators eq.( 4.17| ). 



Beyond the leading order things become more interesting. In particular the independence of 
rapidity does not hold anymore. We will only consider here the double inclusive gluon amplitude 
in detail. The calculation for multi-gluon inclusive is very similar and not more illuminating. At 
order 0(-^-) the double inclusive amplitude has a contribution with a nontrivial dependence on the 
rapidity difference rj — £. This is the new feature of our calculation which is not present neither in 



KLWMIJ nor in JIMWLK limits. It also has not been considered in [2C], as it is formally the next 
to leading order contribution. Similar correlations, however, arise in the BFKL at NLO. 

To calculate the correlated contribution we return to the definition eq.( [4.13 ). The mechanics 



of the calculation is very similar to the single gluon inclusive amplitude. We have to cut two links 
between H's and S''s in eq.( |3.16 ) and eq.( 3.19| ), and settle the cut links with the free transverse 



coordinates and rapidities of the counted gluons. This generates terms of 0(1/ a 2 ), 0{l/a s ) and 
0(1). The former terms add up to the result quoted in eq.( f4.14| ). Terms of 0(1) we neglect and 
concentrate on the terms 0(l/a s ). Some of the terms in the resulting expression do not depend 
on the rapidities of the gluons. We do not have control over this type of terms. The reason is the 
following. As we have discussed in we have not been careful with the subleading terms in Hrpt 
related to relative ordering between b and A. The terms arising from a change of ordering would 
lead to a "virtual correction" to the classical field b of order 5b ~ 0(g). These terms would also 
be present in the expression for Q n and thus would give contribution to our present calculation at 
0(l/a s ). However the structure of these terms is clearly the same as that of the leading 0(l/a^) 
terms, and thus will not depend on rapidity. 

We conclude that at present we do not have control over rapidity independent terms in 0(l/a s ) 
but can unambiguously calculate the rapidity dependent terms to this order. We note that there is 
nothing that fundamentally prevents us from treating the ordering more carefully and thus calculating 
all subleading terms in the nucleus-nucleus scattering case. We choose not to do so in this paper 
since it is a separate question and deserves a careful treatment in its own right. 

We thus concentrate on the terms depending on the rapidity difference r] — £. In principle 

dn(k,p) 



therefore we should be calculating x^^ly, but rather than taking an extra derivative we will 



- 15 - 



simply subtract any term we get which does no depend on rj — £. A straightforward calculation along 
the lines described above gives: 



drjdt; 



\ correlated 



1 

2^ 



ik(x-x)+ip(y-y) 



x([(b- b)N[ s (E^E) {XtTl . m y N ± (b -&) + (&- b)N[.(E^E) {y ^ v) x N ± (b - b) 
+ (b- b)N{-(K^E) (x>v , y>0 y N* ± (b -b) + (b- b)Nl x (tfK) (x>vm y N ± (b - b) 
+ (b- b)N[ s (K^E) (y ^ v) y N* ± (b -b) + (b- b)Nl x (E^K) mx>v) y N ± (b - b) 



)j,S ■ 

(4.20) 

We stress again, that all expressions in this equation and below depend on the single gluon scattering 
matrix S, rather than the matrix R. 

To simplify this expression further we note that the integral of N over rapidity is real, and 
therefore N*(b — b) = N(b — b). Further we note that the matrix E is pure imaginary. Finally we 
have to remember that the operator in eq. (pL20D has to be averaged over the projectile and target 
wave functions. We will assume that the averaging weights are rotationally invariant, and thus the 
result must be invariant under k — > —k; p — > —p, which is equivalent to (x,y — > x,y). Under this 
assumption eq. ((4.20|) becomes 



\correlated 



1 

2^2 



Jk(x-x)+ip(y-y) 



(x,x;y,y) 



drjdc] 

x ( [(6 - b)Nl^E\ x „. u) y N±(b - b) + (6 - b)N^E\ u , x ^ x N ± (b - b) 

+ {b- b)N\_-[E\K - IT) Wo yN±(b - 6) + (6 - b)N[ s [E^K - K*)\ y ^. XyV) y N ± (b - b) 



(4.21) 



)j,s ■ 



We further note that the operators E^E and E*{K — K*) are both Hermitian and real, and therefore 
symmetric. Thus we finally get 



dn(r],k;tj ,p) , 
drjdt] 



{correlated 



7T Z 



Jk(x-x)+ip(y-y) 



(4.22) 



x ( [(6 - b)N{-[E^E] iXjV . m y N ± (b - b) + (b - b)N{-[E\K - #*) We) vN±(b - b)] ) jtS . 
To get a more explicit expression we first note the following two integrals 



1 1 _ 1 § e x ~^ + 1 

\ 1 + Ae^-t 1 + Be x ~Z ~ 2 %e x -i - 1 



hi 



B 
T 



,x- v 



Z 1 + Ae^-t 



1 



+ 



1 



1 - Be x -t + ie 1 - Be x ~t 



le 



+ Ci 



b p \- v 



(4.23) 



|e A -" + 1 



hi 



B 
A 



+ C 2 



The terms C\ and C2 do not depend on the rapidity difference rj — A. They are both formally 
logarithmically divergent when the integration over £ is unrestricted. They do depend in principle 
on r] + A. The dependence on the sum of rapidities however, is simply part of the evolution of this 
observable in the rapidity difference between the projectile and the measured gluons. We will discuss 
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this evolution in the next subsection, but are not interested in it for the current discussion. We thus 
will drop the terms C\ and C2 in the following since they do not depend on the rapidity difference. 
We next use these integrals to calculate the products of matrices E and K that enter eq. ( ff .21 ) : 



[tfE] (XjV . y>x) = { x \L.l±±(X-r,)\y) 



~(x\(l-2l)S(l-2l)(l-2L R ) D , 



n 2 \ 

D 2 e t- J- 



ln 



D 2 L 



n 2 

D 2 



(4.24) 
S\l-2L L )\y) 



n 2 \ 



~(x\(l-2L L )S 



n 1 
D 2 



hi 



n 2 

D 2 



E${K — K*) 



(x,T);y,\) 



(x\(l - 205(1 - 2Z)(1 - 2L fl ) gf e 

%e*-v + 1 



(1 -2^(1-20^(1- 20 \y) 

(4.25) 
S\l-2L L )\y) 



hi 



n 2 

D 2 



7~i 9 - 1 - 



(x\(l-2L L )S 



D 2 



111 



n 2 
D 2 



(l-2L R )(l-2l)SHl-2l)\y). 



There clearly is a nontrivial dependence on the rapidity difference in these expressions. Interestingly 
enough this dependence does not disappear even when the rapidities are far from each other. At 
large values of rapidity difference A — 77 we have 



A)I[A-ij[>1 



E\K — K*) 



l|A-r,|»l 



(4.26) 



{x\2 - (1 - 205(1 - 20(1 - 2L*)St(l - 2L L ) - (1 - 2L L )S{\ - 2L R )(l - 2Z)S+(1 - 2l)\y)\X - rj\ . 

Thus we find that the rapidity correlation between the two gluons in the next to leading order 
does not disappear at large rapidities. Obviously we cannot use these expressions when the rapidity 
difference is too large, since then we have to take into account the rapidity evolution between r\ and 
A. Still r] — A can be taken parametrically of order one, but numerically greater than one. In this 
regime our calculation should be valid and the effect is visible. 

It would be very interesting to understand the physics of the appearance of these correlations 
and their possible implications. It is rather clear that we should understand some of these terms as 
the first correction due to the rapidity evolution between 77 and A, which indeed should be a formally 
subleading correction to the leading result eq.(4.20). At the moment we cannot make any further 
comments on the subject and leave it as a question well worth studying. 

We note however, that not all the terms in eq.( |4.24 ) can be attributed to the evolution. Some 
of them, when integrated over rapidity, contribute directly to the leading order RFT Hamiltonian. 
The Hamiltonian can be represented in terms of the inclusive gluon amplitudes T n in the manner 
analogous to the representation in terms of amplitudes Q n 



Hrft 



J n=V 



■11 



IV. 



_d_ 

dY 



{ki',Vi} 



1 n (ki, ...,k n,T]\, ■■■Tj n )\Y=0 ■ 



(4.27) 
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This picks out the linear in Y piece in the integral of T n , including of course T2, over rapidities. 
Such a contribution does not come from the leading order piece, which is clearly proportional to Y 2 
(for T2) nor from the 7] — £ S> 1 integration region of eq.(4.21), since it is proportional to Y 3 (see 
eq.( 4.26| )). However integrating eq.( [4.21| ) over the two rapidities one will also pick a contribution 
proportional to first power of Y from the region 7/ — ^ ~ 1 where the two gluons are locally correlated 
in rapidity over and above the long range linear rapidity correlation explicit in eq.(4.26). Thus for 
small rapidity differences eq.( 4.21| ) contains local physical effects not related to the rapidity evolution 
between r/ and £. 

Quite generally existence of such short range correlations can be inferred from the structure 
of our calculation of Hrft- Recall that Hrft is obtained by taking the (sum of the) squares of 
exclusive n-gluon production amplitudes Q n , integrating them over the rapidities of all the gluons, 
and picking the piece of the integral that is linear in the length of the total rapidity interval Y. 
Such a linear in Y piece naturally arises from the integral of Q\Qi, since only one rapidity variable 
is integrated over. However for Q n Qn with n > 1 the only way such a linear piece can arise is if 
there is finite excess (or depletion) of probability for all n gluons to be correlated in rapidity within 
a finite interval 771 ~ 772 ~ ... ~ rj n over a totally uncorrelated situation. Since as we have seen earlier 
all QnQn with n > 1 contribute to Hrft, it means that such short range correlations are indeed 
present. Eq.( [4.2l] ) is just a specific example of this correlation. 

We have not attempted to study these local correlations in any detail. It is tempting to speculate, 
however, that they lead to anti-bunching, rather than bunching, thereby depleting the probability 
to emit several gluons close to each other in rapidity. This would then impose sort of a "rapidity 
veto". Such an effect appears in the next to leading order perturbative approaches [28] and has been 
used in the framework of the nonlinear high energy evolution to emulate energy conservation p9|]. 
It would be interesting if the inclusion of Pomeron loop effects discussed here and in |l]] implements 
such a veto automatically. 

We also note that extending our results to include evolution between the rapidities of the observed 
gluons, along the lines of [14] is fairly straightforward and we plan to address this question in near 
future. It would be very interesting to compare this to the approach of the last paper in [20]. 



4.3 Evolution with rapidity 

As we have already stated, we are not going to discuss the evolution with respect to rapidity differ- 
ences between the counted gluons. On the other hand the evolution with respect to the total energy 
is covered by our present derivation in a trivial manner. All the multi-gluon observables discussed 
so far have the form ( |4.1| ) 

0(y,Y) = (0(j,S)) j>s = J [dj}[dS}W y p [j] 0(j,S) W$_ y [S\ (4.28) 

where 0(j,S) = T n and W y is the weight functional for averaging over the projectile degrees of 
freedom, while Wy_ y is the same for target degrees of freedom. The total rapidity of the process is 
Y and the gluons are measured at rapidity y away from the projectile. This expression automatically 
has a factorized structure discussed in f20|| , independently whether the scattering objects are nuclei or 
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dipoles. One can of course ask how does this observable evolve with any of the two rapidity variables 
it depends on: y or Y — y; or in fact with Y at fixed y or Y — y. The answer to this question is 
straightforward given that we know the evolution of the weight function W p derived in the previous 
section. Thus for example the evolution with respect y at fixed Y — y is given by 



d_ 

dy 



0(y,Y-y)\ Y ^ 



[dj] [dS] H RFT [U, R] W p [j] 0(j, S) W?_ y [S] 



(4.29) 



with Hrft given in eq.( |2.34|) . We have not derived directly the evolution of the target weight 
functional. However as was shown in J7j,[27], Lorenz invariance requires the target weight function 
to evolve with the dual Hamiltonian H R ft[Rsi S] where i?g is to S, what R is to U. Namely 



J 



Ta = -- f abc di S^ bd di S dc ; 



R s (x) = V exp{g / dt T 



}. (4.30) 



where j T is the color charge density of the target. Using this we can write the evolution with respect 
to Y at fixed y as 



^0(y,Y-y)\ y 



[dj][dS] W y p [j] 0(j,S) H RFT [R S ,S]W? 



(4.31) 



It was also shown in [g7| that the complete H R ft must be self dual, namely H R ft[U, R] = Hrft[R, U]. 
We have not verified explicitly that Hrft of eq. ( |2.34j ) satisfies the property of self duality. As noted 
in (I]], the technical issue that has to be resolved before we can address this question is the duality 
transformation properties of Ul. This is an interesting question which we plan to address in future. 
Assuming self duality of Hrft, eq.( 4.31| ) can be also written as 



^0(y,Y 



[dj] [dS] W y p [j] 0(j, S) Hrft [S, R s ] W$_ y [S] 



(4.32) 



Combining eqs.( 4.29] ) and ( 4.31 ) one can also write evolution equation with respect to other combi- 
nations of y and Y. Thus if we want to follow the evolution of O with rapidity y keeping the total 
rapidity fixed we have 



^-0(y,Y-y)\ Y 



dj][dS]i H RF T[U,R}W p [j] 0(j,S)W^_ y [S] 



- W p [j]0(j,S) H RFT [S,Rs}W^_ y [S} 



(4.33) 



The factorized structure of this evolution is the same as discussed earlier for dipole-nucleus scattering 
in §§,[14 1 and nucleus-nucleus scattering in |2Qf| . It is rather universal, and one might even say trivial. 
It does not depend on the nature of colliding objects and, by definition, is the generic property of 
any observable that can be represented in the form eq.( 4.28[ ). Of course, not any observable has 
this representation. Examples to the contrary are diffractive observables or observables nonlocal in 
rapidity 

The details of the evolution on the other hand depend on the situation. In the dipole-nucleus 
scattering §§,[14] one does not have to keep the full H R ft in the evolution equation. Instead the 
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Hamiltonian which acts on the projectile weight function can be taken as Hklwmij\3i R]- In the 
nucleus-nucleus scattering, according to p3] the full Hamiltonian reduces to HjjMWLK[U,6/6j]. 



The natural question in this context is, does our formula eq.( 4.33 ) allow to truncate Hrft to the 
JIMWLK form when considering multi-gluon production in nucleus-nucleus scattering. We will 
discuss this question along with some other points in the next section. 

5. Discussion. 

The main result of the present paper is the explicit expression for the single gluon and multi-gluon 
production amplitudes which include Pomeron loop effects. To summarize, in the leading order in 
a s , the n-gluon production amplitude is given by a " semiclassical" expression 

dn(k u kn) = l_f e^^- e 0+-"^(^-*n)([^( Xl ).^(x 1 )]..4 < 4( Xn ). < 4(x n )]) ijS (5.1) 

ar\\...ar\ n tt J (xi,x 1 ;...x n ,x n ) 

where the "classical field" A is an explicit function of the projectile and target fields 

A = N(b-b) = (l-l-L L ) [b L -Sb R ] . (5.2) 
This result reproduces the known expressions for the "dipole- nucleus" scattering. It has the same 



general structure as the expression for nucleus- nucleus scattering of |2C]. The important difference 
between our result and the expressions in |^] is that our expressions for the observables are explicit 
functions of the projectile and target fields, while the procedure of |^] involves further numerical 
solution of dynamical equations of motion to determine the Yang-Mills fields at asymptotically late 
times. In this respect our expression eq.fl5.2j) is the explicit solution of the classical equations of 
motion, needed to apply the formalism of [pCj] . It would indeed be very interesting to check explicitly 
by numerical simulation whether the numerical solution of Yang-Mills equations after the collision 



does reproduce eq.(5.2). We believe this should be the case, once the difference between the gauge 



fixing conditions used in [£0| and in the present paper is accounted for. We note in this connection 
that the eikonal approximation in the symmetric gauge/frame does indeed lead to solution of classical 
equations as the leading perturbative contribution j30|| . Thus we believe that our use of the eikonal 



approximation does not make the accuracy of our approach different from that of [20]. 

We have also shown that two - and higher gluon amplitudes have subleading correlated pieces 
where all gluons at adjacent rapidities are emitted in a correlated way. As we have discussed above, 
the logarithmic part of this correlation is likely to be the first correction due to the rapidity evolution 
between the emitted gluons, while the short range correlation, with the correlation length of order 
unity, is not an evolution effect. These are gluons that upon integration over rapidity contribute di- 
rectly to HftFT- We suspect that these correlations are in fact negative, and thus manifest themselves 
as a sort of rapidity veto due to the Pomeron loop effects. 

Finally we have shown that the multi-gluon observables evolve with rapidity according to Huft 
derived in Q. In this context we again have to discuss the relation of our results to those of ]p0| ], 
Ref . |20f| discusses the evolution of the single and multi gluon production amplitudes with rapidity. 
The rapidity that changes could be either the rapidity of the measured gluons at fixed total energy 
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(evolution with y at fixed Y), or the total rapidity of the process keeping the gluons at mid-rapidity 
(evolution with Y and y so that Y = 2y). In all cases ref. |^] asserts that the observables should 
be evolved with the JIMWLK Hamiltonian. Can we reconcile the two statements? Consider for 
example the evolution 



d dn , 
dy dn 



\Y~y 



[dj}[dS] H RFT [U,R]W y p [j] O g (j,S)W$JS] 



(5.3) 



with the observable O g being the single gluon amplitude O g = AA (|4.2| ) . Since Hrft is a hermitian 
operator, we can also write 



d dn , 
dy dn 



\Y-y 



J [dj][dS] 



W y p [j] H RFT [U,R]O g (j,S) W?_JS] 



(5.4) 



exp{ 5 T<^} in 



Let us now consider the question whether we can expand the operator R(x 
powers of 5/5j when acting on the observable O g . As is clear from eq.(|5.2j), 
function of gb and gb only. The fields gb and gb in turn are functions of gj, as is obvious from the 
classical equations through which they are defined. Thus we have 



the observable O g is a 



An S M 6 
Sj 9 ~ 9 S(gj) 8(gb) 



O 



(5.5) 



Thus each term in the expansion of R when acting on the operator O g brings a factor of a s . The 
same argument goes through for any leading order multi-gluon inclusive amplitude. On the basis of 
this argument we therefore can conclude that for this particular set of observables, one can expand 
the Hamiltonian Hrft to leading order in 5/5 j. As we know, the leading order in this expansion is 



H jiMWLK- Thus, on the face of it, we recover the conclusion of [20] - the evolution hamiltonian for 
these observables can be taken as Hjimwlk- 

The previous line of argument immediately begs two questions. First, how confident are we 
in the robustness of the argument itself. The answer to this question is not completely clear to 
us. Although naively the argument looks correct, there are examples of situations where a similar 
argument fails. Consider the evolution of the dipole-dipole scattering amplitude using the KLWMIJ 
kernel. Ignoring the details irrelevant to the present discussion we can write the evolution equation 
for the S - matrix as 



d 



djda W p [j] K(x, y, z) j, (1 - R) z j y e l £ a ^ W T [a] 



(5.6) 



where /C is the kernel appearing in Hklwmij (@D- Since the target is a dipole, the field a ~ 0(g). 
It then follows by the argument similar to the one given above that when acting on this observable 
( e if x iW^)), the R's in the KLWMIJ Hamiltonian can be expanded to leading order in 5/5 j. The 
leading order expansion is simply the BFKL Hamiltonian. Thus we would conclude that if the target 
is a dipole, we can always use the BFKL Hamiltonian rather than the full KLWMIJ. However we 
know that this is not the case. The BFKL Hamiltonian can be used only for evolution to rapidity of 
order Y ~ — ln^. Further evolution with BFKL Hamiltonian violates unitarity of the amplitude, 
while the full dipole or KLWMIJ Hamiltonian preserves unitarity. Thus even though the formal 
argument about expansion of R at every step of evolution can be made, the cumulative effect of 
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evolution to large enough rapidity is such that the expanded Hamiltonian misses a very important 
physical effect which leads to qualitative change in the evolution. Note that if the projectile is not a 
dipole, but rather a "nucleus" - a state with large dipole density in its wave function, the breakdown 
of Hbfkl happens much earlier, at Y ~ — 1 . In this case we can point out to a distinct physical 
effect that is missed by the expansion - multiple scattering corrections due to scattering of more than 
one dipole of the projectile on the target. Those corrections are important if the projectile contains 
many dipoles, even if the target is dilute. For a dipole projectile multiple scattering corrections 
become important later in the evolution, when the dipole wave function becomes dense, while for 
a nucleus this happens much earlier. In the case of the evolution of multi-gluon amplitudes in the 
nucleus-nucleus scattering discussed above, we do not have similar understanding. Still we think 
that one has to take the argument with a grain of salt. It makes perfect sense to ask how far one 
can evolve these observables in rapidity without encountering a problem of missing some important 
physical effect. We thus would like to advocate caution on this issue. The problem in our view needs 
to be studied further. 

The second question is this. If the argument is indeed correct, does this mean that we can ignore 
the difference between Hjjmwlk and Hrft for all observables, and simply not bother with any of 
the calculations in Q? The answer to this is clearly negative. The fact that the evolution simplifies 
for a certain set of observables, does not mean that it simplifies for all interesting observables. A 
simple example of a similar situation is the dipole-nucleus scattering. In this case the weight function 
of the projectile dipole VF p [j] evolves according to Hklwmij- This evolution is significantly different 
from Hbfkl, since as is well known Hklwmij leads to unitarization of the scattering amplitude, 
while Hbfkl does not. On the other hand, if we consider a single inclusive gluon amplitude, the 
observable itself is quadratic in j (see ||, fljj). Thus when acting on this observable, Hklwmij and 
Hbfkl are identical. The fact that one cannot use Hbfkl to evolve the weight functional is 
another way of saying that there are some interesting observables on which the action of Hk lwmij 
and Hbfkl is not equivalent. As mentioned above, one of such observables is the forward scattering 
amplitude, which is unitarized by the KLWMIJ evolution but not by BFKL. Other examples of such 
observables include various diffractive amplitudes |l3|] and also multi-gluon inclusive amplitudes in 
the case when the rapidity differences between observed gluons are large and the evolution between 
them has to be taken into account [14]. 

We expect the situation to be similar for the case at hand. Clearly, R in Hbft cannot be 
expanded when we calculate the forward scattering amplitude on a nuclear target. In this case the 
observable is exp{ija}. When acting on it, R becomes the matrix S, which is not perturbatively 
close to unity. Thus expansion is not possible. 

We expect that multi-gluon amplitudes with large rapidity differences also do not allow expansion 
of R. The observables associated with these amplitudes are calculable and we hope to present results 
of this calculation soon. However, we can find a hint that all derivatives in R are important, by 
examining the results of the previous section. Consider the correlated term in the double gluon 
inclusive amplitude. We expect that the evolution between the rapidities of the two gluons, when the 



1 We do not claim that the nuclear weight function can be evolved with Hklwmij ■ Physically of course this does 
not make much sense, since the nonlinearities in the wave function are important. We merely point out that if one does 
it as a mathematical exercise, the formal argument about expanding R's breaks down rather quickly. 
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rapidity difference is large, should be given by Hrft, or at least closely linked to it. On the other 
hand, as we have discussed above, the first term in this evolution is likely just the "long range" part 
of the rapidity correlated term eq.( [4.22 ), that is eq.( [4.26| ). Examining eq.( 4.26| ) we see that it is just 
the first order expansion of the denominator of Hrft eq.( 2.34| ) around (1 — 2/)i?^(l — 21). These 
terms would not depend on S if we were to truncate Hrft at Hjimlk, but they clearly depend on 
S in eq.( }4.26| ). This suggests that the factors R cannot be expanded and the complete Hamiltonian 
Hrft is important in the evolution of this observable. We note that the conclusion of ref.[^] is 
different, namely that Hjimwlk is adequate also for evolution of multi-gluon observables with large 
rapidity differences. We feel therefore that the question warrants further study. 

Finally we want to mention that it would be very interesting to understand how to perform the 
averaging over the valence and target fields. This would allow one to calculate physical observables 
like multi-gluon spectra. One can in principle use the McLerran-Venugopalan model [^] to specify 
the projectile and target weight functionals, at least to get a qualitative idea about the behavior 
of the observables. However even then one has to understand how to calculate correlators of three 
(dependent) matrix degrees of freedom Ur, Ul, S, which is far from trivial. This can certainly be 
done numerically, but one would like to be able to understand at least the basics analytically. At the 
moment this is an open question. 
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A. Appendix - Single gluon inclusive amplitude. 

In this appendix we rederive the expression for the single gluon inclusive amplitude without the use 
of the amplitudes Q n , but instead working directly from the definition. 

6 = — (0| & & n a\ a (ri, k) af(ri, k) fit Sfl |0) . (A.l) 

ZTT 

In the parametric range we are interested in, namely when both the colliding objects carry fields of 
order 1/g, we expect the number of produces gluons to be of order l/a s . We will calculate the single 
gluon spectrum only to this order. Let us first of all calculate 

nts+noant^n ( A - 2 ) 

where now the index a stands for all discrete indices as well as momenta. In this equation and the 
following we assume that the rapidity of the gluon operator a(rj, k) is in the infinitesimal rapidity bin 
created by the operator O. Thus in practice we are calculating directly the derivative with respect 
to rapidity. Recall that 

fl = CB. (A.3) 
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Thus 

Qa^ = CB a a B ] C ] = C [Q a pa p + ^ a pa\] C ] . (A.4) 

When acting by the operator C we keep in mind that we are only interested in the leading order 
contribution. Thus we do not need to worry about the action of the operator C on j in O and 
Also recall that 

Ca a C^ sa a a — V2ib a (A. 5) 

In this relation we have neglected terms which are themselves of order one, and are proportional 
to A(x~ = 0). The reason is that we only need to keep terms, which by subsequent application of 
another operator C can be shifted by b, thus generating terms of order 1/g. However A{x~ = 0) 
commutes with C and thus will not generate such contributions. Thus to the required order 

fi a a fit = e a/ 3 ap + $ Q/3 ajj - V2i (8 — $) a p bp[J] . (A.6) 

The subsequent action of the operator S simply rotates all creation/annihilation operators as well as 
all the charge density operators by the single gluon scattering matrix S 

# fi a« fit 5 = 6 a/3 [SJ] Sfr a 7 + <5 Q/3 [SJ] a\ - V2 i N a/3 [SJ] b p [SJ] . (A.7) 

Finally we apply again the transformation with the operator fi. The only relevant part of this 
transformation is the action of the operator C. All the rest, as before does not give a leading order 
contribution 

fit^tfiaafit^fi = VZiQaplSJlSfrhy - V2i^ aP [SJ]Sfs^ - V2iN a(S [SJ}b p [SJ] 

= V2iN aP [SJ] { S Pl 6 7 [ J] - b p [S J] } . ( A.8) 

A similar relation holds for at. Collecting this together we obtain for the single gluon inclusive 
spectrum 

6 g = ^-j e ik(z ~ s) js6[J;x] - 6[5J;x]| N^[SJ;x,z]N[SJ;z,y] jsb[J;y] - 6[5J;y]| (A.9) 
which is the same as derived in the body of the paper (eq. |4.7|) . 
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